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ABSTRACT 

The kinematics, structure, and stellar population properties in the centers of two 
brightest early-type galaxies of the Leo II group, NGC 3607 and NGC 3608, are 
studied by means of integral-field spectroscopy. The kincmatically distinct areas in 
the centers of these galaxies, with radii of 6" and 5" respectively, are found also to be 
chemically distinct. These stellar structures are characterized by enhanced magnesium- 
line strength in the integrated spectra. However, we have not found any mean stel- 
lar age differences between the decoupled cores and their outskirts. Analysis of two- 
dimensional line-of-sight velocity fields reveals systematic turns of kinematical major 
axes near the nuclei of both galaxies; in NGC 3608 the ionized gas rotates in the or- 
thogonal plane with respect to the stellar component rotation. By taking into account 
some morphological features, we conclude that both NGC 3607 and NGC 3608 have 
large triaxial stellar spheroids. We argue that the magnesium-enhanced cores are not 
circumnuclear disks; instead they resemble rather compact triaxial structures which 
may be a cause of formation of polar disks around them - a gaseous one in NGC 3608 
and a stellar-gaseous one in NGC 3607. In the latter galaxy the star formation is 
perhaps still proceeding over the polar disk. 

Key words: Galaxies: individual: NGC 3607 - Galaxies: individual: NGC 3608 
- Galaxies: nuclei - Galaxies: stellar content - Galaxies: kinematics & dynamics - 
Galaxies: evolution 



1 INTRODUCTION 

There exists now a certain contradiction between a homo- 
geneous red-colour appearance of nearby massive early-type 
galaxies implying very ancient epoch of the main star forma- 
tion and prescriptions of the hierarchical concept of galaxy 
formation requiring quite recent merger events. The recent 
merging of disk galaxies possessing significant amount of gas 
must be followed by some secondary, preferably circumnu- 
clear star formation. Galaxy groups represent perhaps the 
best places where consequences of the external-driven evolu- 
tion of early-type galaxies may be found because dense envi- 
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ronments and moderate (with respect to clusters') velocity 
dispersions are favourable for mergers and tidal interactions. 

Recently we have considered Leo I gro up and its mem- 
bers NGC 3379, NGC 3384, and NGC 3368 (jSil'chenko et alj 
2003). A combined analysis of the stellar and gaseous kine- 
matics and of stellar population properties in the centers of 
the galaxies has revealed signatures of synchronous secular 
evolution. In the centers of the galaxies the inner rotation 
axes are aligned, despite the different orientations of the 
outer galactic bodies, and the mean stellar age estimates 
evidence for quasi-simultaneous star formation bursts about 
3 Gyr ago. However, the Leo I group is somewhat unique 
because it po ssesses a superg i ant in tergalactic H I cloud of 
1.7 ■ 10 9 M Q (|Schneiderlll985l ,[T9 8g ). Its shape is a clumpy 
ring with a radius of ~ 100 kpc encircling the galaxy pair 
NGC 3379/NGC 3384, and just the spatial orientation of 
this ring defines the alignment of the inner rotation axes 
of the three galaxies. So it seems probable that the secu- 
lar evolution of the Leo I early-type galaxies is governed by 
the interaction of every galaxy with the intergalactic H I 
ring and not by tidal interactions between the galaxies. In 



2 Afanasiev & Sil'chenko 

Table 1. Global parameters of the galaxies 



NGC 


3607 


3608 


Type (NED 1 ) 


SA(s)0* 


E2 


R 25 , kpc (LEDA 2 ) 


15.2 


10.5 


B° T (RC3 3 ) 


10.79 


11.69 


M S (LEDA) 


-19.84 


-19.46 


(B - V)° (RC3) 
(U - B% (RC3) 


0.92 


0.93 


0.49 


0.40 


V r (NED), km-s" 1 


935 


1253 


Distance, Mpc 




23 4 


Inclination (LEDA) 


34° 


54° 


PA phot (LEDA) 


120° 


75° 


a*, km -s" 1 (LEDA) 


217 


204 



1 NASA/IPAC Extragalactic Database 
2 Lyon-Meudon Extragalactic Database 
3 Third Reference Catalogue of Bright Galaxies 
jTonrv et all J200lf) 



the present paper we consider central early-type galaxies of 
the Leo II group lacking significant masses of neutral hy- 
drogen. It is interesting to look for signatures of the sec- 
ular evolution in the circumnuclear parts of these galax- 
ies. The Leo II group contains 16 ga l axies brighter than 
Bt ~ 16, according to iGiuricin et all |2000t ). From those, 
only 5 are of early type, namely, SO or ellipticals; three 
early- type galaxies, NGC 3605 (compact dwarf elliptical), 
NGC 3607, and NGC 3608 are located in the very center of 
the group within some 60 kpc area. The velocity dispersion 
of the whole group is rathe r moderate, 417 km/s according 
to lMahdavi fc Gellerl l|200ll) . 

The main characteristics of the galaxies to be con- 
sidered in the present paper are given in Table 1. 
NGC 3607 is a giant lenticular galaxy located in the cen- 
ter of the L eo II group. Hot X-ray gas is detected inside 
this galaxy ([Te rashima et al. 2002), as well as around it 
l|Biermann. Kronberg. fc Madorel Il982l ) - all three galax- 
ies in the center of the group are embedde d into the com- 
mon e nvelope of hot gas. An X-ray atlas of iMulchaev et al] 
(2002) demonstrates two maxima of the hot gas concen- 
tration - one at NGC 3607 and another at NGC 3608 - 
which fact forces authors to suggest a recent assembly of 
the Leo II group. Also NGC 3607, though of early type, is 
known to have a rather ext ended ionized-gas disk (IShieldsl 
ll99ll ; lMacchetto et al.H"l996h with a radius of ~ 15" (~ 1.5 
kpc) containi ng a broad dust ring between R = 8'.'4 and 
R = 13'.'2 jSingh et all 1 1994]) . Central structure of the 
neighbouring low-luminosity elliptical galaxy NGC 3608 is 
also very interesting: it has a counterrotating core. From 
their long- slit cross-section taken along the major axis of 
NGC 3608 iJedrzeiewski fc Schechteil (| 19881 ) found that at 
approximately R = 10" (1.1 kpc) the stellar rotation 
changes its sense; the maximum rotation velocity of the in- 
ner core, ac hieved at R ~ 4", w as however rather low, of 
~ 15 km/s. ICarollo et al.l l| 19971 ) inspected high-resolution 
images obtained with the WFPC2/HST through two filters, 
F555W and F814W, for the sample of elliptical galaxies with 
kinematically decoupled cores. In NGC 3608 they have not 
found any evidence for a circumnuclear stellar disk presence. 



Also they claimed that the color V — I is very homogeneous 
over the central parts of all their sample galaxies, the kine- 
matically decoupled cores being indistinctive from the sur- 
rounding stellar bodies. In the particular case of NGC 3608, 
however, we can suspect a break of V — I by less that 0.1 
mag at R ~ 4" from their Fig. Id. Curiously, though either 
ICarollo et af] (|l997h nor another investigators do not find 
any fine morphological substructure in NGC 3608, the global 
properties of the galaxy lo ok unusual: despite its low lum i- 
nosity, it is a slow rotator lljedrzeiewski fc Schechterlll988l ). 
it is boxy (|Bender. Dobereiner. fc Mollenhofll 19881 ). and it 
has a shallow core profile l|Carollo et all 1199% " - a classic 
set of properties of a giant elliptical. General views of both 
galaxies are presented by Figs. [T] and [2] where three various 
fields of view are outlined: the large-scale view of the galax- 
ies, the HST/PC-frame field, and the very central regions 
which we have observed by means of 2D spectroscopy. 

The layout of the paper is as follows. We report our 
observations and other data which we use in Section 2. In 
Section 3 two-dimensional velocity fields obtained by means 
of 2D spectroscopy for the central parts of NGC 3607 and 
NGC 3608 are presented. The radial variations of the stellar 
population properties are analysed in Section 4. Section 5 
gives discussion and our conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

The spectral data which we analyse in this work are obtained 
with two different integral-field spectrographs. Integral-field 
spectroscopy is a rather new approach which was firstly pro- 
posed by Prof. G. Courtes some 20 y ears ago - for a desc rip- 
tion of the instrumental idea see e.g. iBacon et al.l l|l995h . It 
allows to obtain simultaneously a set of spectra in a wide 
spectral range from an extended area on the sky, for exam- 
ple, from a central part of a galaxy. In the spectrographs 
which data we use a 2D array of microlenses provides a 
set of micropupils which are put onto the entry of a spec- 
trograph. After having reduced the full set of spectra we 
obtain a list of the following characteristics for the indi- 
vidual spatial elements: fluxes in continuum and in emis- 
sion lines, line-of-sight velocities, both for stars and ion- 
ized gas, and absorption-line equivalent widths which are 
usuall y expressed as indic es in the well-formulated Lick sys- 
tem l|Worthev e t al. 1994). This list can be transformed into 
two-dimensional maps of the above mentioned characteris- 
tics for the central part of a galaxy which is studied. Be- 
sides the panoramic view benefits, such an approach gives 
an unique opportunity to overlay various 2D distributions 
over each other without any difficulties with positioning. In 
this work we use the data of two 2D spectrographs: the 
fiber-lens Multi-Pupil Fiber Spectrograph (MPFS) at the 
6m telescope of the Special Astrophysical Observatory of the 
Russian Academy of Sciences (SAO RAS) and the interna- 
tional Tiger-mode SAURON at the 4.2m William Herschel 
Telescope at La Palma. SAURON is a private spectrograph, 
but some its raw data are available after expiration for ev- 
erybody from the Isaac Newton Group Archive in the UK 
Astronomy Data Center. 

The last variant of the MPFS became opera- 
tional in the prime focus of the 6m telescope in 1998 



(http:/ /www.sao.ru/hq/lsfvo/devices/mpfs/mpfs.html ). 
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Figure 1. Trje_ j^igtongpjj^, maps for NGC 3607: left - lai^jCgse^Jp yrjajje^ through the filter V (the HautcjP^vcjp^gdata from the 

HYPERLEDA database), intensity is linearly gray-scaled, middle - the HST/PC frame through the filter F547M, the same, right - the 

area corresponding to the field of view of the Multi-Pupil Fiber Spectrograph, with the uncalibrated colour distribution obtained by 

dividing the HST/PC frames through the filters F8f4W and F555W by each other and taking 2.5 logarithms, darker means redder. 
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Figure 2. The ahotpmetric maps for NGC 3608: left - a large-jjcalc red image from DSS, intensity is linearis etay-scaled, middle - 
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the HST/PC frame through the filter F555W, the same right - the area corresponding to the field of view of the Multi-Pupil Fiber 

Spectrograph, with the uncalibrated colour distribution obtained by dividing the HST/PC frames through the filters F8i4W and F555W 

by each other and taking 2.5 logarithms, darker means redder. 



With respect to the previous variant, in the new MPFS 
the field of view has been increased and the common 
spectral range has become larger due to using fibers. In 
2001-2002 the fibers transmitted light from 16 x 15 square 
elements of the galaxy image to the slit of the spectrograph 
together with additional 16 fibers that transmitted the sky 
background light taken 4 arcminutes apart from the galaxy, 
so that the separate sky exposures were not necessary. The 
size of one spatial element was approximately 1" x 1"; a 
CCD TK 1024 x 1024 detector was used before 2003. The 
reciprocal dispersion was 1.35 A per pixel, with a spectral 
resolution of 4-5 A slightly varying over the field of view. 
To calibrate the wavelength scale, we exposed separately 
a spectrum of the hollow cathod lamp filled with helium, 
neon, and argon; an internal accuracy of linearization was 
typically 0.25 A in the green and 0.1 A in the red, and 
additionally we checked the accuracy and an absence of 
systematic velocity shift by measuring strong emission lines 
of the night sky [OI]A5577 and [OI]A6300. We obtained 
the MPFS data in two spectral ranges, the green one, 
4300-5600 A, and the red one, 5900-7200 A. 

The green-band spectra are used to calculate the Lick 
indices H/3, Mgb, Fe5270, and Fe5335 which are suitable to 
determine metallicity, age, and Mg/Fe ratio of old stellar 



populations lw"orthevlll994l ). Following the prescriptions of 
IWorthev! dl994l) we have observed 15 stars from the list 
of IWorthev et all (1 19941 ) during four observational runs to 
calibrate the new MPFS index system onto the standard 
Lick one. After that we have calculated the linear regression 
formulae to transform our index measurements into the Lick 
system. Below we give the regression formulae used by us for 
the index calibration into the standard Lick system: 
H/3 = (0.79 ± 0.09)H/3 M pfs + (0.18 ± 0.13) 
Mgb = (1.04 ± 0.04)Mgb MPJ , s - (0.26 ± 0.16) 
Fe5270 = (0.98 ± 0.05)Fe5270 M PFS - (0.32 ± 0.18) 
Fe5335 = (1.07 ± 0.11)Fe5335 M PFS - (0.78 ± 0.40) 
The comparison of star measurements and the linear calibra- 
tion dependencies are presented in Fig. [3] The rms scatters 
of points near the linear dependencies are less than 0.15 A 
for 3 indices except Fe5335 for which it is 0.27 A, so the in- 
di vidual deviations of s tars are within observational errors 
of IWorthev et all (|l994T ). 

We would like to stress that we do not degrade the 
spectral resolution of our spectra. The Lick indices H/3, Mgb, 
Fe5270, and Fe5335 are absorption- line equivalent widths by 
their definition, and equivalent widths of the spectral lines 
do not depend on spectral resolution: they are integrals, so 
by degrading the spectral resolution you make the lines shal- 
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Figure 3. The correlation of the instrumental absorption-line 
indices obtained with the spectrograph MPFS and the standard 
Lick indices for 15 stars from the list of Worthey et al. (1994). 
The thin straight lines are bisectrices of the quadrants, the dashed 
lines are linear regressions fitted to the dependencies. The mean 
index error for the measurements of the stars by Worthey et al. 
(1994) are indicated for a single star at every plot; our errors are 
lower by an order. 



lower but broader, and integrals remain the same. The only 
restriction is put by the wavelength intervals where the in- 
tegrals are calculated. The Lick index system is arranged to 
include the full line width with the Lick spectral resolution, 
which is in average ~ 8 A, for the majority of galaxies. Both 
the MPFS and the SAURON have the better spectral resolu- 
tion, so in our case the absorption lines under consideration 
are completely within the wavelength intervals prescribed by 
the Lick system, so we need no to make something with our 
spectral resolution. This fact is confirmed by the close coinci- 
dence of the instrumental index measurements of the stan- 
dard stars with the SAUR O N and the tabular index dat a 
from IWorthev et all l|l994h i| Afanasiev fc Sil'chenkol [2002ft . 
However, small variations of the continuum shape caused by 
imperfect spectra calibrations which are individual for every 
spectrograph may give small systematic deviations of the in- 
strumental indices from the Lick index system, and just this 
effect is corrected with the calibration formulae given above. 
In the range of indices used, these corrections do not exceed 
0.4 A for R/3 and Fe5270, and 0.6 A for Fe5335; for Mgb they 
are negligible. For the more detailed description of the Lick 
index measurements with the MP FS and inte r nal a ccuracy 
estimates, see the recent paper by ISil'chenkol (|2006l ) where 
the study of the large sample of lenticular galaxies with the 
MPFS is presented. 

To correct the index measurements for the stellar veloc- 
ity dispersion which is usually substantially non-zero in the 
centers of early-type galaxies, we have fulfilled smoothing of 
the spectrum of the standard star, HD 97907, by a set of run- 
ning Gaussians of various widths; the derived dependencies 



of index corrections on cr* were approximated by polynomi- 
als of 4th order and applied to the measured index values 
before their calibrations into the Lick system. Due to sur- 
face brightness differences between the circumnuclear and 
more outer part of the bulges, the accuracy of the measured 
indices falls from 0.1 A at the center to more than 0.5 A at 
the edges of the multi-pupil frame. To support a constant 
level of the index accuracy along the radii when analysing 
the stellar population properties, we co-add the individual 
spectra within circular rings centered onto a nucleus; the 
typical statistica l error of the az imuthally averaged indices 
is 0.1 A-0.15 A l|Sirchenkoll2006T ). 

The green-range 2D spectroscopic observations are also 
used to cross-correlate galactic elementary spectra with a 
spectrum of a template star, usually of G8III-K3III spec- 
tral type, to obtain in such a way a line-of-sight velocity 
field for the stellar component and a map of stellar velocity 
dispersion. The stars are defocused during the observations 
to take into account the variations of spectral resolution over 
the field of view. For NGC 3607 in particular, we have calcu- 
lated the kinematical maps with three templates of different 
spectral types: HD 73665 (G8III), HD 73710 (G9III), and 
the main component of STF 1947 (K0III). We have checked 
that the kinematical maps obtained for NGC 3607 with three 
different templates do not differ within the errors. The cross- 
correlation peaks are fitted by Gaussians; the benefits of this 
approach as opposed to the more popular Fourier Correla- 
tion Q uotient (FCQ) are thoroughly discussed by iBottemal 
(1988); we would like to stress that the cross-correlation 
method is not so sensitive to template mismatch as FCQ 
or direct-fitting methods. The red spectral range contains 
the strongest optical emission line for LINERS, [NIL1A6583, 
so it is used to derive line-of-sight velocity fields for the ion- 
ized gas by calculating emission-line baricenter positions for 
[NII]A6583. The accuracy of elementary velocity measure- 
ments, both for stars and ionized gas, is about 10 km/s. 

The second integral-field spectrograph which data we 
use in this work is a new instrument, the SAURON, oper- 
ated at the 4.2m William Herschel Tele scope (WHT ) on L a 
Palma - for its detailed description see l Bacon et al.l (|200lf ). 
We have taken the data for NGC 3608 from the open ING 
Archive of the UK Astronomy Data Centre. Briefly, the field 
of view of this instrument is 41" x 33" with the spatial el- 
ement size of / .'94 x 0"94. The sky background taken 2 ar- 
cminutes from the center of the galaxy is exposed simulta- 
neously with the target. The fixed spectral range is 4800- 
5400 A, the reciprocal dispersion is 1.11 A-1.21 A varying 
from the left to the right edge of the frame, and the spec- 
tral resolution is about of 4 A. The comparison spectrum is 
neon one, and the linearization is made by a polynomial of 
the 2nd order with an accuracy of 0.07 A. The index sys- 
tem is checked by using stars from the list of lWorthev et all 
(1994) which have been observed during the same observa- 
tional run. The regressions fitting the index system calibra- 
tion of the February- 1999 run wh en NGC 3608 has been 
obser ved are shown in our paper (|Afanasiev fc Sil'chenkol 
2002). The relations between instrumental and standard- 
system indices are very close to the y — x relation so no 
corrections are needed to calibrate them into the standard 
Lick system. The stellar velocity dispersion effect has been 
corrected in the same manner as for the MPFS data. While 
to prepare the azimuthally averaged index profiles from the 
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MPFS data we co-added the spectra in the rings, to pre- 
pare the SAURON azimuthally-averaged index data which 
have higher signal-to-noise ratios, we averaged the measured 
individual-element indices over the same rings so the at- 
tached error bars are the errors of the means formally cal- 
culated - they are all below 0.03 A. 

The full list of the exposures made for NGC 3607 and 
NGC 3608 with two 2D spectrographs is given in Table 2. 
'BTA' means 'Bolshoi Telescope Azimuthal'ny', it is an of- 
ficial name of the 6m telescope. 

To analyse the structure of the galaxies and to refine a 
kinematical analysis, for both galaxies we have retrieved the 
WFPC2/HST image data from the HST Archive. NGC 3607 
was observed in the frame of the program of A. Phillips (ID 
5999) on nuclei of SOs, and NGC 3608 was observed in the 
frame of the program of M. Franx on kinematically decou- 
pled nuclei (ID 5454). For NGC 3607 we have also used 
the large-scale images taken from the database HYPER- 
CAT/FITS Archive (PI Ph. Prugniel). The details of the 
photometric observations are given in Table 3. 

All the data, spectral and photometric, except the data 
obtained with the MPFS, have been reduced with the soft- 
ware produced by Dr. V.V. Vlasyuk in the Special Astro- 
physical Observatory l|Vlasvuklll993l ). Primary reduction of 
the data obtained with the MPFS was done within IDL with 
a software created by one of us (V.L.A.). The Lick indices 
were calculated with our own FORTRAN program as well as 
by using the FORTRAN program of Dr. A. Vazdekis which 
provides also the calculation of statistical errors for the in- 
dices. 



3 STELLAR AND GASEOUS KINEMATICS 
AND THE STRUCTURE OF THE CENTRAL 
KILOPARSECS 

Since the integral-field spectroscopy provides us with two- 
dimensional line-of-sight velocity fields, we are able now to 
analyse both character of rotation and central structure of 
the galaxies. If we have an axisymmetric mass distribution 
and mean rotation on circular orbits around the symmetry 
axis, the direction of maximum central line-of-sight velocity 
gradient (we shall call it 'kinematical major axis') should 
coincide with the line of nodes of the main symmetry 
plane as well as the photometric major axis should do. In 
the case of a triaxial shape of potential the isovelocities 
would tend to align w i th th e principal axis of the ellipsoid 
dVauterin fc Deionghd [19971 simulations of bars and refer- 
ences therein). In a general case of triaxial potential the 
kinematical and photometric major axes diverge, showing 
turns with respect to the line of nodes in opposite senses if 
the main axis of the tri axial potential is not strictly aligned 
with the line of nodes (|Monnet. Bacon, fc Emsellemll 19921 ; 
iMoiseev fc Mustsevovl [ioOOh . In the simplest case of thin- 
disk rotation we have a convenient analytical expression 
for the azimuthal dependence of the central line-of-sight 
velocity gradients within the area of solid-body rotation: 

dvr/dr = usinicos(PA — PAq), 

where ui is the deprojected central angular rotation velocity, 
i is the inclination of the rotation plane, and PAo is the ori- 
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Figure 4. The azimuthal dependencies of the central stellar 
line-of-sight velocity gradients taken in several radial bins for 
NGC 3607 according to the MPFS two-dimensional velocity mea- 
surements. 



entation of the line of nodes of the rotation plane. In more 
complex three-dimensional cases the validity of the formula 
above is not so evident; however the observed dependen- 
cies of dv r /dr on PA also look harmonic as one can see in 
Figs. 2] and Perhaps, it results from symmetry properties 
that vertical and radial velocity components vanish after 
integrating along the line of sight. At least it was shown 
that axisymmetric oblate spheroids demonstrate ze ro line- 
of-sig ht velocity gradients along their minor axes (|Binnevl 
1985) and th e maximum rotation along their major axes 
(|Baconll 19851 ). So by fitting azimuthal variations of the cen- 
tral line-of-sight velocity gradients by a cosine law, we can 
check the axisymmetry of the mass (stellar?) distribution 
and to determine the orientation of the kinematical major 
axis as its phase and the central angular rotation velocity as 
its amplitude. It is our main tool for kinematical analysis. 



3.1 NGC 3607 

Figure [6] presents mostly kinematical maps for the central 
part of NGC 3607. Both stars and ionized gas demonstrate 
regular, almost rigid-body rotation with their kinematical 
major axes close to the photometric major axis. The gas ro- 
tation is twice faster than that of the stars: the central angu- 
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Table 2. 2D spectroscopy of the galaxies studied 



Date 


Galaxy 


Exposure 


Configuration 


Field 


Spectral range 


Seeing (FWHM*) 


29 Apr 01 


NGC 3607 


45 min 


BTA/MPFS+CCD 1024 x 1024 


16" 


x 15" 


4200-5600 A 


2'.' 3 


09 Mar 02 


NGC 3607 


45 min 


BTA/MPFS+CCD 1024 x 1024 


16" 


x 15" 


5800-7200 A 


2'.' 8 


09 Mar 02 


NGC 3608 


60 min 


BTA/MPFS+CCD 1024 x 1024 


16" 


x 15" 


5800-7200 A 


2'.' 8 


20 Feb 99 


NGC 3608 


120 min 


WHT/SAURON+CCD 2k x 4fc 


33" 


x 41" 


4800-5400 A 


1"4 



Table 3. Photometric data on NGC 3607 and NGC 3608 



Date 


Galaxy 


Telescope 


Filter 


Exposure 


Seeing 


Scale 


7 May 1994 


NGC 


3608 


WFPC2/HST 


F555W 


500 s 


O'.'l 


0'.'045 


7 May 1994 


NGC 


3608 


WFPC2/HST 


F814W 


230 s 


O'.'l 


0'.'045 


6 Nov 1994 


NGC 


3607 


WFPC2/HST 


F547M 


260 s 


O'.'l 


0'.'045 


27 Jan 1996 


NGC 


3607 


WFPC2/HST 


F555W 


160 s 


O'.'l 


0'.'045 


27 Jan 1996 


NGC 


3607 


WFPC2/HST 


F814W 


160 s 


O'.'l 


0'.'045 


6 Feb 1999 


NGC 


3607 


1.2-OHP 


B 


300 s 


6'.'3 


0'.'686 


6 Feb 1999 


NGC 


3607 


1.2-OHP 


V 


120 s 


6'.'0 


0'.'686 


6 Feb 1999 


NGC 


3607 


1.2-OHP 


R 


120 s 


6'.'0 


tf'686 


6 Feb 1999 


NGC 


3607 


1.2-OHP 


I 


120 s 


5'.' 7 


a'686 



20 - 
10 - 
- 
-10 - 
-20 - 
20 - 

1 o - 



-10 - 
-20 - 
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Figure 5. The azimuthal dependencies of the central stellar 
line-of-sight velocity gradients taken in several radial bins for 
NGC 3608 according to the SAURON two-dimensional velocity 
measurements. 



lar rotation velocity is u> ga s sini ~ 50 km/s/arcsec whereas 
Tatars sin i « 28 km/s/arcsec. Let us to note that if the incli- 
nation given by LEDA for the symmetry galactic plane and 
derived from the large-scale disk axis ratio, 34° (Table 1) , is 
correct, it means that the central gas rotates at a speed of 
800 km/s/kpe. However, below we should recognize that it 
is rather difficult to determine precise orientation(s) of the 
rotation plane(s) in the center of NGC 3607. The stellar ve- 
locity dispersion has a peak near the nucleus. In our data 
this stellar velocity dispersion maximum area is not exactly 
centered onto the nucleus that may be a real effect due to 
the asymmetric dust distribution in the center of the galaxy; 
this maximum area seems to be extended in the direction of 
the isophote minor axis. 

To check our 2D mapping of the line-of-sight velocities, 
we have simulated one-dimensional cross-sections of the ve- 
locity fields of Fig. [5] a band aligned with the major axis 
has been overposed onto the maps and the velocity and 
dispersion velocity values inside this band have been de- 
rived and averaged over the 1" x 2" elements. We have com- 
pared the simulated LOS velocity and stellar velocity disper- 
sion profiles with the long-slit data along the major axis of 
NGC 3607 which we have found in the literature (Fig.[7]and 
[8]). One can see that the agreement of the velocity profiles is 
rather good. However, the literature data do not support the 
shift of the velocity dispersion peak with respect to the pho- 
tometric center. Also, perhaps the slope of the ionized gas 
velocity profile in our data may be slightly underestimated 
due to moderate seeing quality of the observations of March 
2002. The maximum gas rotation velocity, ~ 250 km/s at 
R — 5" — 10" which may be transformed into 450 km/s under 
i = 34°, seems to be enormous; but this result is confirmed 
by recent CO observations of NGC 3607 according to which 
2.3 ■ 10 8 Mq of the m olecular gas within R = 5" rotates at 
a speed of 250 km/s l|Welch fc Sagell2002i ). Interestingly, in 
the integrated molecular-line profile only a receding horn is 
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Figure 6. The line-of-sight velocity fields of the stellar component (top left) and of the ionized gas according to the measurements of 
[NIIJA6583 (bottom left), and the stellar velocity dispersion map (top right) and the [Nil] emission line intensity distribution in arbitrary 
units (bottom right) in the central part of NGC 3607. The continuum, green (A5000 A) for the stars and red (A6400 A) for the gas, is 
presented by isophotes. 



well seen - the asymmetry is quite similar to that of the 
Fig. [6] (the lower left). The stellar velocity profile of Fig. [7] 
reveals a prominent rotation velocity maximum at R ~ 6" 
after which the rotation velocity drops by a factor of ~ 3 
toward R = 20" — 25". Such behaviour may be treated as a 
kinematically distinct core in the center of NGC 3607; let us 
to note that the maximum rotation radius coincides exactly 
with the border of the magnesium-enhanced structure (see 
the next Section). 

The enormous rotation velocity, which we have mea- 
sured in the center of NGC 3607 and which implies a mass 
concentration of more than 10 Af@ within R ~ 0.5 kpc, 
provokes some doubts about a circular (axisymmetric) char- 
acter of rotation. We can check it by comparing kinematical 
and photometric major axis orientations; it is done in Fig. [9] 
Our first impression from inspecting Fig. [6] is that the kine- 



matical and photometric major axes coincide; the Fig. [9] 
confirms this impression by quantifying it. At R > 3", or 
more exactly in the radius range of R = 3" — 6" , within the 
decoupled core area, where the rotation is still almost solid- 
body, the kinematical major axes of the ionized gas and 
of the stars, as well as the photometric major axes of the 
continuum brightness distribution are all close to the global 
line of nodes of the galaxy, PA = 300° (120°). This fact 
agrees with an axisymmetric character of rotation with the 
spin orthogonal to the main symmetry plane of the galaxy 
- a classic configuration corresponding to the SAO-type of 
NGC 3607. But inside R = 2" ALL THE AXES - both the 
kinematical ones and the photometric one - seem to turn to- 
ward larger PAs, the turn of the gas kinematical major axis 
being the most prominent. We have estimated an orientation 
of the innermost isophotes of Ha+[NII] emission distribution 
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Figure 7. The comparison of the line-of-sight velocity profiles 
simulated along the major axis by using our 2D velocity fields for 
the stars and ionized gas in NGC 3607 with the literature data 
of Fisher (1997) and of Caon et al. (2000). The slit width used in 
the simulations is 2". 
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Figure 8. The comparison of the stellar velocity dispersion pro- 
files simulated along the major axis by using our 2D velocity fields 
for the stars and ionized gas in NGC 3607 with the literature data 
of Fisher (1997) and of Caon et al. (2000). The slit width used in 
the simulations is 2". 
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Figure 9. Isophote major axis position angle compared to the 
orientations of the kinematical major axes (see the text) for the 
stars and ionized gas in the center of NGC 3607. The errors of the 
kinematical major axes determination are estimated as 1—1.5 deg. 
The line of nodes determined from the outermost disk isophote 
orientation, PA = 300°, is traced by a dashed line. 
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Figure 10. Isophote characteristics in the full radius range of 
NGC 3607. 



in Fig. If of lShieldsl (|l99l[ ) - it is PA w 325°, very close to 
our measurements of the [Nil] emission-line brightness dis- 
tribution elongation at R > 3". The gray-scaled em ission 
brightness map provided by Macchctt o et al.l |l990)) con- 
firms also the turn of the innermost emission isophotes to- 
ward larger PAs. The concordance of the photometric major 
axis provided by the HST data and the kinematical major 
axis of the ionized gas near PA w 325° -330° at R w 1" -2" 
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implies axisymmetry. But the rotation plane of the inner- 
most, R < 2", part of the galaxy, and so that of circum- 
nuclear disk probably, appears to be inclined to the main 
symmetry plane of NGC 3607, while the more outer gas ro- 
tates in the global galactic plane together with the stars. The 
latter fact rejects an external origin of the decoupled central 
spin; there must be own intrinsic peculiarities of NGC 3607 
to provide such kinematical appearance. Figure \W\ shows ra- 
dial variations of the isophote characteristics in NGC 3607 
at larger scales than Fig. [9] does. The isophote major axis 
position angle remains almost constant all along the radius. 
But the ellipticity behaviour is indeed strange: it reaches a 
main maximum of 0.3 at the border of the decoupled core, 
at R — 6", then it reaches a secondary maximum of 0.2 at 
R » 20", and after that it falls to 1 - b/a « 0.1 at R « 60". 
The formal inclination of the global galactic disk, i = 34° 
(LEDA), corresponds only to 1 — b/a = 0.17 by assuming in- 
finitely thin disk. Let us also to note that the radius of the 
minimal ellipticity, R m 60", coincides with the boundary 
betwe en the bulge and the global disk of NG C 3607: accord- 
ing to lBaggett, Baggett, fc Anderson! (|l998l 1. in NGC 3607 
the bulge dominates at R < 60" and the exponential disk 
- at larger radii. The ellipticity behaviour at R > 20" is 
roughly consistent in all filters, so it is not an effect of dust. 
We must conclude that the bulge of NGC 3607 cannot be 
an axisymmetric (oblate) spheroid if to require its symmetry 
plane to coincide with the disk plane, otherwise its visible 
ellipticity should be smaller than that of the more thin disk. 
It has to be triaxial at a scale of a few kiloparsec. 
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Figure 11. The line-of-sight velocity fields of the stellar com- 
ponent (top left) and of the ionized gas according to the mea- 
surements of [NIIJA6583 (bottom left), and the stellar velocity 
dispersion map (top right) and the [Nil] emission line intensity 
distribution (bottom right) in the central part of NGC 3608. The 
continuum, green (A5000 A) for the stars and red (A6400 A) for 
the gas, is presented by isophotes. 



3.2 NGC 3608 

To analyse the stellar kinematics of the central part of 
NGC 3608 we have used the SAURON data reduced by us 
with our tools. Recently the stella r kinematical m a ps for 
NGC 3608 have been published by lEmsellem et al.l (|2004l ^ 
together with the maps of LOSVD higher moments /13 and 
/14; the SAURON team has applied its own method to an- 
alyze these data which is quite sophisticated. The quanti- 
tative agreement of the LOS velocity fields obtained by us 
and by lEmsellem et all (|2004l ) is perfect; our stellar veloc- 
ity dispersion measurements being qualitatively similar are 
systematically lower by about of 40 km/s. In the further dis- 
cussion we use mainly the morphology of the stellar velocity 
dispersion distribution, so the precise absolute values of a, 
are not so important in this particular work. 

In Fig. [11] one can see the comparison of the LOS ve- 
locity fields for the stars and for the ionized gas. It reveals 
indeed a central counterrotating stellar subsystem with the 
kinematical major axis close to the major axis of the inner 
iso photes - it is the kinemat i cally decoupled core reported 
bv lJedrzeiewski fc Schechterl (1988). The ionized-gas veloc- 
ity field (Fig. 1111 bottom) reveals a strong twist of isove- 
locites over the whole measured velocity field. The ampli- 
tude of the circumnuclear azimuthal variations of the gas 
line-of-sight velocity gradient is 73 ± 15 km/s/arcsec - even 
higher than in NGC 3607. The phase of the cosine curve fit- 
ting the azimuthal variations of the central gas line-of-sight 
velocity gradients within R = 2" is shifted by ~ 70° with 
respect to the line of nodes of the stellar rotation. As the 
nucleus of NGC 3608 is inactive and radio-quiet, we would 
disregard a hypothesis of outflow and would claim rather 
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Figure 12. Isophote major axis position angles compared to the 
orientations of the kinematical major axes (see the text) for the 
stars and ionized gas in the center of NGC 3608. The line of nodes 
determined from the outermost isophote orientation, PA = 255° , 
is traced by a dashed line. The data for the stars are zoomed 
separately at the top plot to demonstrate the kinematic major- 
axis turn. The errors of the kinematical major axes determination 
are estimated as 1—1.5 deg. 
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a kind of strongly inclined gaseous ring with a radius of 
~ 200 pc. In Fig. [I2]we compare the orientations of the kine- 
matical and photometric major axes, for the stars and for 
the ionized-gas emission [NHJA6583. The kinematical and 
photometric major axes of the stellar component coincide 
that implies an axisymmetric rotation. However both ma- 
jor axes begin to deviate from the PA of outer isophotes, 
PAo — 255°, in a systematic way when approaching the 
very center. Meantime, the major axes of the [Nil] emission 
brightness distribution and of the ionized-gas velocity field 
are turned strongly with respect to the stellar configuration 
at 7? < 4", or within the kinematically decoupled area. The 
spatial resolution of our red-band MPFS exposures is rather 
moderate, so the measurements inside R — 2'.' 5 are strongly 
affected by the seeing effect. From the published HST data 
analysis l|Tomita et al.ll2000t ) we know that in the center of 
NGC 3608 the dust ring with the radius of 0'.'5 is roughly 
aligned with the outer-isophote orientation, PAdust ~ 80° 
(or 260°) (see also our Fig. [2] right plot). Since the dust and 
the gas are related, the orientation of the circumnuclear gas 
plane must also be close to the line of nodes, so one may ex- 
trapolate our measurements back to the line of nodes inside 
R ~ 1". Strong twist of both the photometric and kinemati- 
cal major axes for the gas component at R = 1" — 4", though 
being obviously smoothed by the seeing effect, implies nev- 
erthless an existence of a kind of polar gas rotation in the 
above-mentioned radius range. 

Figures [13] and [14] demonstrate simulated major- and 
minor-axis cross-sections of the stellar LOS velocity and ve- 
locity dispersion fields in the center of NGC 3608 in com- 
parison with published long-slit data. The agreement of our 
cuts with the literature long-slit data is satisfactory though 
the systematic shift of our stellar velocity dispersion mea- 
surements by ~ —30 40 km/s can be also seen. We would 

note that systematic differences of a few dozen kilometre per 
second are still common for stellar velocity dispersion mea- 
surements. For example, in t he rec ent study of 48 early-type 
galaxies bv lEmsellem et all (|2004f ) the individual differences 
between the measured r e /10-aperture stellar velocity disper- 
sions and the literature data reach 40 km/s (see their Fig. 3). 
Partly this uncertainty may be due to using different meth- 
ods to derive kinematical information from the absorption- 
line spectra: various methods are sensitive to signal-to-noise 
ratio and to template mismatch by various degree. However 
?) have shown that even using the same method over the 
data of similar quality does not save the matter: the stellar 
velocity dispersion estimates made over different spectral 
ranges may differ by 20-30 km/s. So we would not refer to 
the exact absolute values of <r», but rather would analyse 
their relative variations over the field of view. 

Besides the overall check of the quality of our velocity 
measurements, Figs. [T3l and 1141 confirm once more that the 
edge of the kinematically decoupled stellar core is located 
at R = 4" and, as we shall see in the next Section, coin- 
cides roughly with the border of the magnesium-decoupled 
core. But the kinematical analysis does not leave a pos- 
sibility to treat the decoupled core as a compact stellar 
disk. The low v/cr ratio, ~ 0.1, and a pattern of the stel- 
lar velocity dispersion distribution - an extended maximum 
clearly seen in the stellar velocity dispersion map (Fig. 111! 
top right) - allow to suggest a triaxiality of the decoupled 
core. In such configuration the polar rotation of the ionized 
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Figure 13. The comparison of the line-of-sight velocity profiles 
simulated along the major and minor axes by using the SAURON 
2D LOS velocity field for the stars in the center of NGC 3608 with 
the literature data of Jedrzejewski and Schechter (1988). The slit 
width used in the simulations is 2". 



gas becomes explicable: we often meet circumnuclear polar 
gaseous disks in gal axies with triaxial bulges or within bars 
- e. g. in NGC 2841 (tSiT'chenko. Vlasvuk. fc Burenkovlll997l: 
Afanasiev fe Sil'chenkol Il999h. in NGC 6340 (ISil'chenkd 
20001 ), in NGC 7280 d Afanasiev fc Sil'chenkd l200Ch . in 
NGC 4548 l|Sil'chenkdl2002l l. etc. 



4 CHEMICALLY DECOUPLED CORES IN 
NGC 3607 AND NGC 3608 

4.1 NGC 3607 

Figure [15] presents Lick index maps for the central part of 
NGC 3607 constructed from the MPFS 2D spectral data. 
The Mgb map reveals a certain presence of chemically de- 
coupled core; this finding is not very unexpected because 
NGC 3607 has been earlier listed by us as a g ood candi- 
date for possessing a chemically distinct nucleus (Sil'chenko 
Il994h . Interestingly, the (Fe) =(Fe5270+Fe5335)/2 and H/3 
absorption line index maps look rather homogeneous and do 
not demonstrate any sharp features. Perhaps, one can note 
a very shallow unresolved (Fe) minimum near the photo- 
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Figure 15. The MPFS index maps for NGC 3607; < Fe >=(Fe5270+Fe5335)/2. The green (A5000 A) continuum is overlaid by isophotes 
in three coloured plots. At the right top the Mgb index distribution smoothed strongly with the 2D Gaussian of FWHM = 3'.' 5 is plotted 
by isolines to show an orientation of the chemically decoupled structure; here the green continuum is gray-scaled 
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Figure 14. The comparison of the stellar velocity dispersion 
profiles simulated along the major and minor axes by using the 
SAURON 2D velocity dispersion field for the stars in the center of 
NGC 3608 with the literature data of Jedrzejewski and Schechter 
(1988), Halliday et al. (2001), and Proctor (2002) The slit width 
used in the simulations is 2". 



metric center of the galaxy and similarly shallow extended 
H/3 minimum aligned with the major axis of the isophotes; 
the latter feature may be associated to the circumnuclear 
ionized-gas disk mentioned in the Introduction and is prob- 
ably resulted from the emission contamination of the H/3 
absorption line. Interestingly, all the central peculiar fea- 
tures of the index distributions are centered slightly be- 
low and to the left from the photometric center showing 
the same shift as the stellar velocity dispersion. However, 
both (Fe) =(Fe5270+Fe5335)/2 and H/3 absorption- line in- 
dex minima are marginal, and the Mgb peak is much more 
prominent than the other details. It seems to be resolved; 
moreover, it seems to be elongated along the minor axis of 
the isophotes. To derive the orientation of the magnesium- 
enhanced core, we have smoothed strongly the Mgb map 
and have drawn isolines of the smoothed Mgb distribution 
(Fig. 1151 top right). Indeed, the magnesium-enhanced core is 
elongated in PA ~ 46° though the Mgb isolines do not look 
as quite symmetric ovals. The orientation of the magnesium- 
enhanced area hints that we see some stellar substructure 
with its alignment along the minor axis of the isophotes 
- some kind of minibar or so. But in fact, the Mgb-index 
distribution morphology alone is not enough to make a con- 
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Figure 16. The azimuthally-averaged index profiles for 
NGC 3607 in comparison with the literature data. The minor- 
axis long-slit cross-section by Proctor (2002) is plotted. The long- 
slit data taken along the major axis by Fisher et al. (1996) are 
also compared to the major-axis index profiles simulated from 
our 2D index maps with the digital slit width of 2". The central 
2arcsec X 4" aperture measurements by Trager et al. (1998) are 
plotted as connected squares at R = 0" — 2". 



elusion about the nature of the magnesium-enhanced core. 
The visible elongation of the magnesium-enhanced core may 
be a result of a magnesium depression along the major axis, 
similar to that of the H/3 index, which may be caused by 
a very young stellar disk aligned with the major axis. To 
make a definite conclusion, a further analysis of the stellar 
population age in the area of magnesium-enhanced core is 
needed. 

To quantify peculiarities of the absorption-line in- 
dex distributions in the center of NGC 3607, we 
have calculated azimuthally averaged radial profiles of 
the indices and have compared them to the 'zero- 
dimensional' an d one-d i mensi o nal spectral data pub- 
lished earlier by IProctorl d2002h. fTrager et al.1 (| 199^1 . and 
iFisher. Franx. fc Illingworthl (|l996l ) (Fig. [HI). Though the 
long-slit data were taken along the major or minor axis of 
NGC 3607 and so were not obliged to agree exactly with our 
azimuthally averaged measurements, they confirm the over- 
all shape of the index radial dependencies: a rather flat iron 
profile, an H/3 profile with a shallow minimum in the center, 
and a magnesium-distinct core with the radius of 5" — 6". A 
moderate, by 0.2 A-0.4 A, systematic difference of our H/3 
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Figure 17. 'Index-index' < Fe > vs Mgb diagnostic diagram for 
the azimuthally averaged Lick indices in the center of NGC 3607 
taken along the radius with the step of 1" (open circles connected 
by a thin line, the nucleus being marked as 'nuc'). The errors of 
the azimuthally averaged MPFS indices are about 0.1 A. The 
models of Thomas et al. (2003) for [Mg/Fe]=0.0, +0.3, and +0.5 
are plotted as a reference frame; the small signs connected by 
pointed, dashed, and point-dashed lines represent stellar popu- 
lation models of equal ages; the metallicities for the models are 
+0.67, +0.35, 0.0, -0.33, -1.35, and -2.25 if one takes the signs 
from top to bottom 
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Figure 18. 'Index-index' H/3 vs [MgFe], [MgFe] = (Mgb < Fe > 
) 1//2 , age-diagnostic diagram for the azimuthally averaged Lick 
indices in the center of NGC 3607 taken along the radius with 
the step of 1" (open circles connected by a dot-dashed line, the 
nucleus being marked by 'nuc'). The errors of the azimuthally 
averaged MPFS indices are about 0.1 A. The models of Thomas 
et al. (2003) for [Mg/Fe]=0.0 and +0.3 are plotted as a reference 
frame; the small signs connected by solid and dashed lines repre- 
sent stellar population models of equal ages; the metallicities for 
the models are +0.67, +0.35, 0.0, -0.33, -1.35, and -2.25, if one 
takes the signs from the right to the left. 



and Mgb profiles with the data of iFisher et all (|l996t ) can 
also be noted. To check if these shifts may be due to a lack 
of azimuthal symmetry, we have simulated one-dimensional 
major- axis cross-sections of our index maps and have also 
plotted the simulated results in Fig. 1161 Indeed, the simu- 
lated major-axis Mgb profile has a much better agreement 
with the Fisher's et al. data than the azimuthally averaged 
measurements. As for the H/3 index profiles, the simulated 
major- axis profile and the azimuthally averaged one have 
coincided almost perfectly, implying an orientation of the 
ionized-gas disk far from edge-on. We must conclude that 
the systematic shift of ^0. 2 A between o ur measurements 
of the H/3 index and those of lProctorl (|2002t ) and lFisher et al] 
( 1996) exists really. Meant i me, th e centered aperture mea- 
surements by iTrager et all (|l998l ) which are a Lick-system 
etalon confirm rather our calibration. 

Due to progress in evolutionary population synthesis 
during last years, we can now estimate mean stellar pop- 
ulation characteristics by confronting different absorption- 
line indices to each other. Figure [17] presents an 'iron-vs- 
magnesium' diagram where we compare our azimuthally av- 
eraged M PFS data for the central part of NGC 36 07 with the 
models of iThomas. Maraston. fe Bender! ll 20031 ) which are 
calculated for several values of magnesium-to-iron ratio. The 
loci of the models of various [Mg/Fe] are well separated on 
the diagram, so from inspecting Fig.[l7]we can conclude that 
the magnesium-to-iron ratio in the center of NGC 3607 is 
certainly above the solar one. But whereas in the circumnu- 
clear region the [Mg/Fe] is between zero and +0.3, perhaps, 
closer to +0.2, the unresolved nucleus is outstanding with 
its [Mg/Fe] « +0.4. We must note that, as Fig. \W\ demon- 
strates, the previous Lick index measurements for NGC 3607 



do not confirm this jump of [Mg/Fe] at R — 0" and imply 
rather constant iron-index behaviour along the radius. Since 
our estimates of the nuclear indices are confined to a single 
spatial element, we cannot insist that the low Fe5270 value 
in the unresolved nucleus is not a random artifact. So over 
all the center of NGC 3607 we assume [Mg/Fe] to be between 
+0.2 and +0.4. 

Several inves t igator s, e.g. iGonzalej (|l993h or 
iTerlevich" fe Forbes! (|2002h . noted that in order to overcome 
an effect of the non-solar magnesium-to-iron ratio when 
determining a mean age of stellar population, one must 
use a combined index [MgFe]= (Mgb(Fe)) 1/2 . So in Fig. [TH] 
we confront [MgFe] to H/3 and compare our azimuthally 
averaged data for NGC 3607 with two sets of models - 
for [Mg/Fe] = and for [Mg/Fe] = +0.3. As it has been 
promised by our choice of the metal-line index, the age 
estimates are robust to the varying choice of the model 
[Mg/Fe]; they are equal to about 10-12 Gyr. The nucleus 
is again outstanding, probably, due to the Balmer emission 
concentration in the very center of the galaxy and obvious 
H/3 emission contamination. But beyond the nucleus the age 
gradient along the radius is non-detectable: at the diagram 
'H/3, [MgFe]' the index variations occur along the model 
sequence of equal age, namely, of T — 10 ± 2 Gyr. We may 
suspect that the real mean age of the stellar population in 
NGC 3607 is slightly lower than 10-12 Gyr, beca use the 
ioniz e d-gas disk is known to e xtend up to R « 15" (I Shields! 
Il99ll : iMacchetto et al] |l996) and therefore the Balmer 
emission may contaminate the absorption-line H/3 index 
in some degree all the way. However, our measurements 
of the [OIII]A5007 emission- line equivalent width in the 
nucleus and in the summed spectrum of its nearest outskirts 



14 Afanasiev & Sil'chenko 



(R = 4" - 7") have given 0.7 A for the former and A 
for the latter, so the effect of contamination of the H/3 
index by emission must be strong only for the nucleus. We 
try to correct it here by using the measurements of the 
Hq-em ission equivalent width bv lHo. Filippenko. fc Sargent! 
(1 19971 ). To correct the measured Lick index H/3 for the 
emission we use the well-known fact that the emission line 
Ha is always much stronger than the H/3 and the absorp- 
tion line Ha is always weaker than the H/3 one, so the 
equivalent width of the emission line Ha can be measured 
more precisely than that of H/3. iHo et all (|l997h obtained 
EW (Hat emis) = 2.03 A for NGC 3607 by subtracting 
pure-absorption template from the aperture nuclear spec- 
tra. Let us to note here that in the nucleus of NGC 3608, 
another our target galaxy, they saw only a marginal Ha 
emission, less that EW = 0.4 A, and no emission lines in 
the green spectral range. The well-established and minimum 
possible intensity ratio Ha/H/3 is known for the radiative 
excitation of gas by OB-stars ('HII-region'-type excitation); 
it is 2 .5-2.7. For sho ck excitation of the gas it is much 
larger. IHo et all ([19971 ) have classified the nuclear emission 
spectrum of NGC 3607 a s a L INER, not a 'HII-type' 
nucleus. Stasinska & Sodrc ( 20011 ) have analysed a large 
heterogeneous sample of integrated spectra of galaxies 
of various morphological types and have found a good 
correlation £VU(H/3 emis) — 0. 25 E W (Ha emis); we use 
just this relation to calculate EW (H/3 emis) , which is in 
fact the correction of the Lick index H/3 for the emission. 
After it has been corrected for the emission contamination 
of the H/3 index, the nucleus of NGC 3607 has settled to the 
age sequence of 12 Gyr in the Fig. 1181 and any age difference 
between the nucleus and its outskirts has disappeared. We 
would like to stress that any significant age break is also 
absent at the border between the magnesium-enhanced core 
of NGC 3607 and its 'bulge', at R = 5" - 6". 

4.2 NGC 3608 

Figure [T5] presents the Lick index maps for NGC 3608, sim- 
ilar to those presented for NGC 3607 in Fig. 1151 but calcu- 
lated from the SAURON data, so instead of the combined 
iron index (Fe) =(Fe5270+Fe5335)/2 we give here only the 
Fe5270 surface distribution. We can note a qualitative re- 
semblance of the index distributions in both galaxies: just 
as we have seen in NGC 3607, in its elliptical neighbour the 
magnesium index demonstrates a prominent extended peak 
in the center, and the iron index shows a flat, homogeneous 
distribution. However, the isolines of the smoothed Mgb dis- 
tribution (Fig. 1191 top right) are elongated in PA « 75°, 
or along the global major axis of NGC 3608 (see Table 1), 
conversely to NGC 3607 where the magnesium isolines are 
aligned with the minor axis of the isophotes. The ellipticity 
of the Mgb isocontours within R m 5" is very high, 0.33, 
which value exceeds the isophote ellipticity anywhere in the 
galaxy, 1 — b/a < 0.2. Recently the SAURON team has pub- 
lished their own results of the reduction of the same data 
l|Kuntschner et ail [2006 ). Though not quantified in detail, 
these results have also implied the extended magnesium- 
enhanced region in the center of NGC 3608 with the el- 
lipticity of Mgb-isolines higher than the ellipticity of the 
isophotes. 

In Fig. [20] we present azimuthally averaged Lick in- 
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Figure 20. The SAURON azimuthally averaged Lick index pro- 
files for NGC 3608 in comparison with the literature data. The 
minor-axis long-slit data of Proctor (2002) are plotted. The aper- 
ture measurements by Trager et al.(1998) and by Kuntschner et 
al. (2001) are also plotted as connected squares and triangles at 
K=0"-2". 



dex profiles and their comparison to the literature data. 
The shapes of the profiles are qualitatively similar to those 
in NGC 3607 (Fig. I16p . Perhaps, the magnesium-enhanced 
core is somewhat more compact: we would put its border at 
R « 5". The break of the Mgb-profile slope at this radius is 
clearly seen: between R = 0" and R = 4'.' 7 the Mgb index 
changes by 0.8 A, whereas beyond R = 5" approximation of 
the azimuthally-averaged measurements by a linear fit gives 
the slope of —0.011 ± 0.010 A per arcsec, or almost negli- 
gible, and an intercept of 4.27 ± 0.008 A - compare to the 
measured central Mgb value of 5.09 A. Let us to remind 
that the maximum rotation veloc ity of the counterrotating 
core is also achieved near R « 4" ijjedrzeiewski fc Schechterl 
1988), so in the case of NGC 3608 we deal with a dynam- 
ically and chemically decoupled co re, just as in the case o f 
another elliptical galaxy NGC 4365 (jBender fc Surmall992l ). 
The Fe5270 measurements are in g o od ag reement with the 



Ihe rteozYU measurements are m good agreement with the 
minor-axis long-s lit data of Proctorl ( 20021 ) and with central 
aperture data of iTrager et af] (| 19981 )7 so we may conclude 
that in the SAURON data of February 1999 for NGC 3608 
the systematic offset of the Fe5270 cali bration by 0.4 A seen 
in the later data e.g. for NGC 3384 i|de Zeeuw et al.ll2002l : 
ISil'chenko et al.ll2003l ). is probably absent. However, there 
may be a small systematic shift, by some +0.2 A, of the H/3 
calibration. 

With these azimuthally averaged Lick index profiles, we 
would try to estimate radial variations of the mean charac- 
teristics of stellar population in the center of NGC 3608 
by using index-index diagrams. In Fig. [21] (top) we com- 
pare the SAURON index radial variations with the mod- 
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Figure 19. The SAURON Lick index maps for NGC 3608. The green (A5000 A) continuum is overlaid by isophotes in three coloured 
plots. At the right top the Mgb index distribution smoothed strongly with the 2D Gaussian of FWHM = 3'.' 5 is plotted by isolines to 
show an orientation of the chemically decoupled structure; here the green continuum is gray-scaled 



els of lThomas et alj l| 20031 ) for [Mg/Fe]= +0.3 at the dia- 
gram 'Fe5270, Mgb'. The agreement looks rather good: both 
chemically decoupled core and its outskirts have [Mg/Fe]~ 

+0.2 h0.3. So to determine a mean stellar age in this 

galaxy, we must take the models with just this magnesium 
overabundance. Figure ^T\(bottom) shows the comparison of 
the azimuthall y averaged Lick index data for NGC 3608 with 
the models of lThomas et all l| 2003h for [Mg/Fe]= +0.3 at 
the diagram 'H/3, Mgb'. Again, as in the case of NGC 3607, 
we do not detect any age gradient along the radius: the data 
sequence follows the model sequence of equal age, this time 
it lies between the model age sequences of T = 5 Gyr and 



T — 8 Gyr. Since the emission lines are very weak in the 
center of NGC 3608 (in particular, we have failed to detect 
noticeable [OIII]A5007 emission line in the green SAURON 
spectra), we think that this age estimate, T = 5 — 7 Gyr, is 
close to the real value. Similarly to NGC 3607, NGC 3608 
does not reveal any age break at the transition point from 
the magnesium-enhanced core to the rest of the galaxy. The 
moderate metallicity difference of AlogZ ft; +0.2 is how- 
ever clearly detected between the nucleus, R = 0", and the 
'non-core' stellar body, R > 4" (Fig. [21] bottom) . 

Let us to note that if the systematic shift of the H/3 
calibration by some 0.2 A is real, this would result in the 
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Figure 21. 'Index-index' diagnostic diagrams for the azimuthally 
averaged Lick indices in the center of NGC 3608 taken along 
the radius with the step of CK'94 (solid points connected by a 
dashed line, the nucleus being marked by 'nuc'); (top) Fe5270 
vs Mgb diagram, with the models of Thomas et al. (2003) for 
[Mg/Fe]= +0.3, and (bottom) H/3 vs Mgb, with the same model 
set. Small signs in the top plot connected by solid lines belong to 
stellar population sequences of equal ages, particularly, of 5 Gyr, 
8 Gyr, 12 Gyr, and 15 Gyr from top to bottom; their metallicities 
are +0.35 and 0.00, if one takes the signs from the right to the 
left 



age un derestimation by some 4 Gyr. For example, IProctorl 
( 2002) has found for the center of NGC 3608 the mean stel- 
lar age of ~ 9 Gyr, although the abundance characteristics, 
[Mg/Fe]« +0.2 - +0.3 and [m/H]= +0.4 in the nucleus and 
+0.2 in the 'bulge', given by him coincide with our esti- 
m ates. Similarly , if for NG C 3607 the H/3 index calibration 
by iFisher et~afl (|l996T ) and IProctorl ((2002) is more correct 
than ours, it would mean that we overestimate the mean 
stellar age in the center of NGC 3607 by the same 4 Gyr. So 
we do not insist on the absolute values of the age estimates 
obtained here. The main result of this consideration is that 
the relative ages do not show noticeable variations along 
the radii within some 10" from the nuclei in both galaxies, 
and particularly, that there are no age drops or rises at the 
borders of the magnesium-enhanced cores. 



5 DISCUSSION 

After a thorough study of the central regions of the bright- 
est galaxies of the Leo II group, the nearest to the center of 



the group, NGC 3607 and NGC 3608, we note a kinemati- 
cally distinct area in the center of the former and confirm a 
presence of the counterrotating core in the latter. The pres- 
ence of the coun terrotating stellar core in NGC 3608 had 
been reported by IJedrzeiewski fc Schechterl (|l988h : a fall of 
the rotation velocity after the circum nuclear maxim um at 
R = 6" in NGC 3607 was noted by iFisherl <|l997l ). Both 
kinematically decoupled cores have appeared to be also dis- 
tinguished by a higher magnesium absorption-line strength, 
but do not differ from their outskirts as concerning Fe5270 
((Fe)) or H/3. 

The most popular hypothesis of the kinematically de- 
coupled core or igin is a hypothesis of smaller elliptical 
satell i te sinking dKormendvll 19841 ; IJedrzeiewski fc Schechterl 
1 19881 ; iBalcells fc Quinnl 1990T I: the smaller ellipticals have 
denser cores that can survive during accretion onto a gi- 
ant galaxy. Later, when lrBender fc Surma! l|l992f ) have found 
that the kinematically decoupled core of NGC 4365 is also 
magnesium-enhanced, the locations of the decoupled rota- 
tion and the magnesium enhancement being coincident, they 
have suggested that the decoupled core of NGC 4365 has 
been formed in a secondary star formation burst after a dis- 
sipative minor merger and subsequent gas accumulation in 
a circumnuclear disk. Taking in mind that NGC 3607 and 
NGC 3608 resemble NGC 4365 as concerning kinematically 
decoupled, magnesium-enhanced cores, for our two galax- 
ies we would however reject both hypotheses on the decou- 
pled core origin mentioned above. The decoupled cores of 
NGC 3607 and NGC 3608 cannot be accreted small ellip- 
ticals because the galaxies have peaks of the stellar veloc- 
ity dispersion and of the magnesium-line strength just in 
the centers of the decoupled cores. Meatime in the frame 
of the dissipationless accretion scenario the satellite light 
must dominate in the center of the merger product and 
so must demonstrate a lower metallicity in the center than 
that of the giant hosts, and the ve locity dispersion must 
have a dip fealcells fc Quinnl 1 1990T ). But the scenario of 
iBender fc Surmal l|l992Tl ~ is also inapplicable: in the frame of 
the dissipative merger hypothesis, the formation of metal- 
enriched stars must proceed in a circumnuclear disk, and 
in the cases of NGC 3607 and NGC 3608 we have proved 
that the magnesium-enhanced structures are not disks. In 
both galaxies the areas of the maximum stellar velocity 
dispersion are elongated. Dynamical simulations show that 
in a triaxial potential the high velocity dispersion areas 
must have an oval shape and to b e aligned with the model 
bars (jVauterin fc Deionghd fl9 97f). The association of the 
magnesium-enhanced areas in NGC 3607 and NGC 3608 
with some ellipsoidal bodies is confirmed also by absence of 
disk-like rotation, when tracing the photometric major axes 
of these structures. In NGC 3607 the magnesium-enhanced 
circumnuclear feature is elongated directly to the kinemat- 
ical major axes of the stars and of the ionized gas - evi- 
dently, it is orthogonal to the circumnuclear disk major axis. 
In NGC 3608 the stellar rotation, though proceeding with 
the expected spin orientation, is too slow for the visible el- 
lipticity of the magnesium-index isolines. If we attribute the 
Mgb-isoline ellipticity, 0.33, to a hypothetical circumnuclear 
disk, then the minimum v max /a ratio appropriate to oblate 
spheroids or disks would be 0.7 (llllingworthlll977h , and we 
have ~ 0.1 in the center of the galaxy. The ionized gas ro- 
tates in the polar plane with respect to the stellar rotation 
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in the very center of the galaxy. So we conclude that the 
chemically decoupled cores in NGC 3607 and NGC 3608 
look like small magnesium-enriched triaxial structures with 
some kinds of polar disks around them - only gaseous one in 
NGC 3608 and stellar-gaseous one in NGC 3607. We think 
that such a configuration may be possible if occurs inside a 
global triaxial spheroid. 

Thoug h NGC 3607 is classified as SO , some investi- 
gators, e.g. iJedrzeiewski fe SchechteH (|l988T l. thought it to 
be an elliptical capturing its diffuse outer envelope from 
NGC 3608 during their interaction. Perhaps, there exists 
a principal difference between formation of the decoupled 
cores in elliptical and disk galaxies. In disk galaxies the 
chemically decoupled structures are indeed circumnuclear 
stellar disks in the most cases. When we have found such 
struc tures in the lent icular galaxies NGC 1 023 (Sil'chcnko 
1 19991 ) and NGC 3384 jSil'chenko et al.ll2003l ). we have noted 
that they are distinguished by a higher strength of both mag- 
nesium and iron absorption lines, the iron-index enhance- 
ment being spatially more extended than that of the Mgb. 
We then concluded that the secondary star formation bursts 
having produced the decoupled cores were brief and effective 
in the unresolved nuclei, resulting in their increased [Mg/Fe] , 
and prolonged during at least a few Gyr in the circumnu- 
clear disks. In the centers of NGC 3607 and NGC 3608 which 
may accrete gas from a 3D quasi-spherical hot halo, and not 
from an outer fast-rotating disk nor from a supergiant flat 
gaseous ring as it was the case of Leo I group, there were no 
perhaps initial favourable conditions for disk formation. Un- 
der spheroidal burst geometry, its duration was short, and 
only magnesium enhancement had a time to be imprinted 
into the secondary stellar populations. 

An interesting detail is the absence of the mean age 
difference between the decoupled cores and the main bod- 
ies of NGC 3607 and NGC 3608. In disk galaxies, including 
the above mentioned lenticulars NGC 1023 and NGC 3384, 
the chemically decoupled cores are always 'younger' than 
the bulges, the mean (luminosity-averaged) stellar age of 
the nuclei being 3-7 Gyr versus 10-15 Gyr of the bulges. It 
implies that the secondary star formation bursts producing 
the decoupled cores in the disk galaxies are rather recent. 
In NGC 3607 and NGC 3608 the decoupled cores do not 
seem to be disti nguished by a y ounger mean age of the stel- 
lar population. iProctorl (|2002h who studied radial profiles 
of the stellar population properties along the minor axes of 
NGC 3607 and NGC 3608 with the long-slit spectra has 
not found a n y age breaks at R = 5" — 6" too. Finally, 
iDavies et al.l 1 200 ll) who study 2D distributions of the stel- 
lar population properties in the center of NGC 4365, the el- 
liptical galaxy with kinematically decoupled core, have not 
also found any age difference between the core and its broad 
outskirts, both being rather old, T > 12 Gyr. Such age 'ho- 
mogeneity' is in some contradiction with the hypothesis of 
the secondary nuclear star formation burst having produced 
the decoupled cores. However it can be explained by an older 
age of the decoupled cores in elliptical galaxies: as the pho- 
tometric and spectral evolution becomes very slow at larger 
ages, one needs a very high precision of the measured Lick 
indices to detect an age difference of a few Gyr when both 
decoupled cores and their outskirts are older than, say, 8 
Gyr. When a larger statistics on the ages of decoupled cores 



in elliptical galaxies become available, this idea would be 
checked more carefully. 
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